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Abstract 

We present a novel approach called the intermediate rotating wave approximation (IRWA), which employs a 
time-averaging method to encapsulate the dynamics of light-matter interaction from strong to ultrastrong coupling 
regime. In contrast to the ordinary rotating wave approximation, this method addresses the co-rotating and 
counter-rotating terms separately to trace their physical consequences individually, and thus establishes the continuity 
between the Jaynes-Cummings model and the quantum Rabi model. We investigate IRWA in near resonance and large 
detuning cases. Our IRWA not only agrees well with both models in their respective coupling strengths, but also offers 
a good explanation for their differences. 

Keywords: Ultrastrong coupling regime, Jaynes-Cummings model, Rabi model. Intermediate rotating wave 
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1. Introduction 

The quantum Rabi model (QRM), which describes 
the interaction between a qubit and a quantized 
harmonic oscillator (bosonic mode) HI, is written as 

Hgabi = fro a O~z + tur) r a ' « + %CT x (a + a 1 ), (1) 

where a (a f ) represents the bosonic annihilation 
(creation) operator of the electromagnetic field mode, 
oj r is the corresponding frequency; to a is the transition 
frequency of the qubit, cr ( - (i = x,y,z) are the 
corresponding Pauli operators; and g is the dipole 
interaction strength. The QRM has been widely applied 
in modern physics, ranging from condensed-matter 
physics ED, atomic physics S3 to quantum optics 
(4), such as cavity QED [0 and circuit QED (3 |7| 
systems. Given its great importance, the QRM has 
been studied extensively using various methods m 
|9)- Despite all those studies, the exact solution of 
QRM was only obtained by Braak recently ifTOl . This 
analytical solution, however, is in the form of composite 
transcendental function defined in power series. The 
search for simpler analytical solution of generalized 

Preprint submitted to Elsevier 


QRM with more physical insights thus continued ini 
HD- For example, the Bogoliubov-type transformation 
M M IH is used to diagonalize the Hamiltonian to 
gain solutions and properties of the model 03 ED- 
The quantum Rabi model can be further simplified 
into the renowned Jaynes-Cummings model (JCM) 
m by rotating wave approximation (RWA) provided 
that the coupling strength is sufficiently weak (g <s 
mm{u> r , u> a }), and the detuning is small enough (\oj a - 
io r \ <s u> r + u> a ). In the interaction picture, the rotating 
terms cr+a and <r a ’ oscillate slowly with phase factor 
of exp [±i(oj a -a> r )f], whereas the two “counter-rotating” 
terms cr+a 1 ' and cr_a oscillate rapidly with phase factor 
of exp[+i(cu,. + oj a )t]. Together with weak coupling 
condition, one can separate the time scales and discard 
the fast-oscillating terms mm p.354], thus obtaining 
the Jaynes-Cummings Hamiltonian 

Hjc = -ho> a cr z + hu> r a^a + hg(cr+a + cr_a*), (2) 

which has simple analytical solutions. 

Enhancement and tunability of light-matter 
interaction is crucial not only for fundamental studies 
of cavity/circuit QED but also for their applications 
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in quantum information processing. Three different 
coupling regimes can be defined based on the basic 
frequency scales of the system. In the weak coupling 
regime (g <$; {y, /<•}, k and y being the loss rate of the 
photon and the emitters’ excitation), the discrete density 
of photonic states modifies the radiative lifetime of the 
quantum emitters (Purcell effect) (20). Strong coupling 
regime is achieved when {j,k} <k g <s min {u) r ,a> a } 
EB p.432], such that quantum emitters absorb and 
spontaneously re-emit a photon many times before 
dissipation becomes effective. This strong coupling 
regime has been investigated in various systems, 
ranging from atoms 0, through quantum dots (QD) 
ll22l to Cooper-pair boxes |6j]. In these conventional 
QED experiments, the system is operating in either 
weak coupling regime or strong coupling regime. 
Therefore, the RWA, which leads to JCM from 
QRM, is very well justified, and the JCM captures a 
wealth of physical phenomena in conventional QED 
systems comprehensively. With recent advances of 
new technologies, the ultrastrong coupling regime has 
become experimentally accessible in semiconductor 
l23l |241 l25ll and superconducting systems f7] [26). 
In this so-called ultrastrong coupling regime E3, 
the coupling strength becomes comparable to the 
frequency of the resonator, g/a> r > 0.1. Therefore, 
the routinely invoked RWA and the JCM break down, 
and the systems dynamics become governed by the 
QRM. This novel unexplored physics has opened up 
new research interest in applications of the QRM. 
Since then, considerable progress has been made and 
fascinating phenomenon have been predicted, such as 
photon blockade (28l . nonclassical state generation 
||29) . breakdown of the standard master equation [13011. 
and ultrafast two-qubit quantum gate operations ED- 
One interesting observation is that the co-rotating 
terms and counter-rotating terms in the QRM affect 
the system in a different manner depending on the 
coupling regime. In light of this, one could gain better 
physical intuition of the continuity between the JCM 
and the QRM by treating the coupling strength of the 
co-rotating terms and counter-rotating terms separately 
Il32l [331. In this paper, we seek to understand the 
emergence of the counter-rotating terms from the JCM 
to the QRM by resorting to the time-averaging method 
0 P- 353], which also helps us to keep track of the 
time scale involved in the dynamics. With this method, 
we deploy a form of approximation, which we term 
as intermediate RWA (IRWA). The basic idea of this 
approximation is that, instead of going to the limit of 
either RWA or non-RWA, we use the time-averaged 
coupling strength in the interaction Hamiltonian. We 


present the general formalism and apply the IRWA into 
two specific situations: the near resonance case and the 
dispersive (large detuning) case. 

The paper is organized as follows. In Sec. [2] we first 
give a brief review of the time-averaging approach and 
introduce the IRWA. In Sec. [3j the IRWA is used to 
study the energy levels of the system by perturbation 
theory in the near resonance case with increasing 
coupling strength. In Sec. [4] the dynamics of the system 
is investigated in the dispersive case with IRWA for both 
single- and multi-qubit case. We summarize our results 
in Sec. 0 

2. Time-averaging and intermediate RWA 

2. 1. Time-averaging function 

The slow and fast time scales in a dynamical system 
can be separated explicitly by means of a temporal 
filtering operation. The time average of a function is 
defined by the convolution 

X oo /-*oo 

= dt'zcr(t')f(t+f), (3) 

OO sj — OO 

where the averaging function m{t) is positive, m{t) > 0, 
even, tzr(f) = m(-t), and normalized, f dt ui(t) = 1 

0 p.353]. The weighting function m(t) has a temporal 

r r°° ^jl/2 

width t = I dfz 3 r(f)f“ < 00, which washes out 
oscillation with period smaller than r. 

A simple example of such function is a Gaussian 
function 

1 

m(t) = —— e . (4) 

r \2n 

It is more convenient to work in the domain 
of frequency, by using the Fourier transformed 
time-averaging function, 

X OO 

dt m(f)e iL> ', (5) 

00 

which is real and even, K(-lj) — K(io) = K*(u>), and 
has a finite width of u>k ~ 1/t. The Fourier transform 
of the convolution Eq.([3]l is just the product of the 
individual Fourier transforms: 


F(oj) = K(aj)F(a>). (6) 

K{u>) is also called the cut-off function, which acts on 
F(co) in such a way that F(u>) is essentially unchanged 
for small frequencies, u> <K u>k, whereas frequencies 
larger than the width cj » u> K are strongly suppressed. 
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2.2. Time-averaged Hamiltonian in intermediate RWA 

In order to apply the time-averaging function to the 
QRM, we impose the condition 

g<zm K (7) 

where the cut-off frequency <ok is chosen in such a way 
that the state interaction-picture state | i/'t/)) is essentially 
constant over the averaging interval, i.e. \if/{t)) ~ \i//(t)). 
Upon time-averaging, the Schrodinger equation of the 
QRM in the interaction picture can be written as |4] p. 
354] 

Q 

ih— 1 ^( 0 ) = Hint |, ( 8 ) 

ot 

where time-averaged Hamiltonian H in ,(t) reads 

= H intr (t) + H lntar (t), (9) 

Hi„t, r (t) = frgr (ao-+e lAt + a^cr_e~ ,A '') , (10) 

) = tigar [acr-e -1 ^ + aV+e' 2 ') . (11) 

Here, the time-averaged coupling strengths for 

co-rotating term g r and for counter-rotating term 
g ar , are modified by the cut-off functions, such that 
g r = K(±A)g and g ar = K(±L)g , with A = to a - co r 
and X = to a + to r . This guarantees that the co-rotating 
terms cr+a, cr_a 1 ' and the counter-rotating terms cr + a\ 
cr-a contribute differently to the dynamics of the 
system, depending on the separation of the frequency 
scales. Notice that since the cut-off frequency io K is 
coupling strength dependent (c.f. Eq.(j7]i), the cut-off 
function K(to) thus is a function of both g and to, 
i.e. K(co, co K (g))- 

Going back to Schrodinger picture, we then have the 
time-averaged quantum Rabi model as 

H Ra bi = H 0 + H r + H ar , (12) 

where 

H r = Hg r X+, with X ± = acr + ± a + cr_; (13) 

H ar = frgarY+, with Y ± = acr_ + a 1 cr+. (14) 

The condition of |A| < X is generally satisfied in 
cavity/circuit QED systems, but this does not give a 
justification for us to neglect the contribution of the 
counter-rotating Hamiltonian H ar , and the QRM is still 
needed to describe the system. In the case of RWA, the 
sufficiently weak coupling condition, g min {to r ,to a }, 
allows us to separate the frequency scales by 

g <s to K <s min{co r ,u) a } < X. (15) 



Figure 1: The ratio of the time-averaged coupling strengths g ar for 
counter-rotating terms to g r for co-rotating terms as a function of 
coupling strength gloj r , with Gaussian weighting function of width 
= 10g for small detuning, A = 0.0 \(i) r . As the coupling strength 
is further increased, even though the near resonance condition is 
satisfied, the weak coupling condition is no longer respected and 
hence, RWA is not applicable any longer. The inset is a zoom of the 
region of g/co r between 0 and 0.1. 

Therefore, when the coupling strength is weak 
compared to the free energy of the system, 
counter-rotating contribution is negligible since the 
cut-off function K(Z) is vanishingly small. In Fig. [T] we 
show the ratio between the two time-averaged coupling 
strengths g ar and g r as a function of the normalized 
coupling strength g/to r , with Gaussian weighting 
function of width a> K = 10# for small qubit-resonator 
detuning A = 0.0 \to r . We note that, with small coupling 
strength, i.e. g/co r < 0.05, the time-averaged coupling 
strength g ar for counter-rotating terms is negligible 
for small detuning, which agrees well with the RWA 
conditions. We then arrive at the time-averaged JCM in 
Schrodinger picture, 

Hjc = + hto r a' i a + hg r X + . (16) 

However, as the coupling strength is getting larger, the 
contribution of the counter-rotating terms is increasing 
and hence, needs to be included to describe the 
dynamics correctly. 

3. Near resonance case in intermediate RWA 

In the near resonance case of |A| <s min{ai a , to r ] <s 
oj fl + (o r , the time-averaged coupling strength for 
counter-rotating terms is much smaller than the 
time-averaged coupling strength for co-rotating 
terms, g ar <s g r . Thus, we take the time-averaged 
counter-rotating Hamiltonian in Eq.([T4|) as a 
perturbation to the time-averaged Jaynes-Cummings 
Hamiltonian in Eq.([T6| and apply the non-degenerate 
stationary perturbation theory to obtain the 
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second-order modification for energy ll34l p.249]. 
The results can be written as 


£n = £ (0) + £ (1 ) + E (2) 

°’g 0,g 0,g ^ c 0,g’ 

E„, ± = E ( „°l + E ( " ± + E ( n % (n= 1,2,3...), 


(17) 

(18) 


where 


£l 0) = -Ehu) a , 

o ,g o 
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Figure 2: (Color online) The system spectrum as a function of 
the coupling strength g/oj r with the time-averaged counter-rotating 
Hamiltonian of the IRWA in second perturbation theory (black solid 
lines) compared to the JCM (blue dashed-dotted lines) and the QRM 
(red dashed lines) for A = 0 and o>k = 10 g. 


with C„ = cos 0 n , S n = sin#„, and 0„ = 

arctan( 2 g r V« + 1 /A). 

The eigenenergies of the ground and low-lying 
excited states as a function of coupling strength g / u> r are 
plotted in Fig. [2] which are obtained by our approach of 
IRWA (black solid lines), the JCM (blue dashed-dotted 
lines) and the QRM (red dashed lines). As shown in 
figure Fig. [2] the JCM curves deviate from the QRM 
curves for increasing coupling strength. The curves 
of the second-order perturbation theory in IRWA agree 
with the numerical results of the QRM for ultrastrong 
coupling regime of g/u> r up to about 0.3. 

4. Dispersive case in intermediate RWA 

In this section, we will study the dynamics of the 
Rabi Hamiltonian in the dispersive limit, where the 
qubit and the resonator are far detuned compared 
to the coupling strength g <s |A|. The dispersive 
regime is of practical interests with applications in 
many cavity/circuit QED systems, such as quantum 
non-demolition measurement of the qubit lf35l . parity 
measurement of the two- or multi-qubit state ll36l . 
and quantum gate operations E2. However, most 
of these applications were studied under the condition 
of strong coupling regime, where RWA still holds. 
In 2009, Zueco et al. generalized the studies of 
dispersive Hamiltonian to the ultrastrong coupling 
regime ED- Although the effect of counter-rotating 
terms is merely quantitative in single qubit case, the 
treatment beyond RWA gives rise to a qualitatively 
different effective model for multi-qubit scenario. Using 


4 













































































the time-averaging functions in Sec. [2] to keep track 
of both the co-rotating and counter-rotating terms in 
the Hamiltonian, we can gain a better insight of their 
contributions and hence, the transition from the JCM to 
the QRM. 

4.1. Dispersive Regime in One Qubit 

In the dispersive limit, where the coupling strength is 
much lesser than the qubit-resonator detuning, 

g « |A|, (29) 


number n = (a’a), is the ac-Stark shift j 6 j and g 2 /A is 
the Lamb shift mm. On the other hand, given the 
fact that the resulting Hamiltonian H d r commutes with 
cr z , i.e. [H d , cr,] = 0, it allows quantum non-demolition 
measurement since the qubit’s state will not be changed 
upon the evolution of the system. Hence, the state 
of the qubit can be inferred by probing the resonator 
frequency. 

In ultrastrong coupling regime, where either or both 
of the RWA conditions are violated, we have the 
following inequalites instead. 


the time-averaged QRM in Eq.([T2| can be transformed 
to 


Hi =UH Rahl U f 
h 

~-m a a- z + 


1 1 
+ ggrgar I T + ^ 


Hu r a f a + — + ^ 7 -j [cr. (la^a + l)J 

. (ci^~ + a 2 j 


up to second order in A = g r / A and A = g flr /X via the 
unitary transformation 


U = exp[TX_ + AK_], (31) 

X _ and Y are defined in Eq.(p~3|> and Eq.fB}. In strong 
coupling regime with RWA, we have the following 
inequalities 


g «: |A| «: u) K «: X, (32) 

which encapsulate the dispersive limit and the RWA 
conditions (near resonance and weak coupling limit) 
while respecting the time-averaging condition Eq.(j7]). 
When these inequalities are satisfied such that g ar - 
0 , g r = g. the counter-rotating terms can be safely 
discarded. This gives rise to 

H d r = r ^[w a + ^ j cr- + h ^u r + ^-cr z j o' a, (33) 

where the oscillator frequency is shifted as 

to,. -> to,. ,. = (jo r ± g 2 /A, (34) 

depending on the state of the qubit. Similarly, the level 
separation of the qubit is shifted to 

,(g 2 ,-g 2 t \ 

co a -» w a ,r = oo a + I — +2— a'a\, (35) 

depends on the number of photons in the resonator. 
The term la^ag 2 /A, which is linear in the mean photon 


g <s |A| < X «: u K , 


(36) 


and all terms in Eq.(30i will be retained, and the 


effective Hamiltonian then reads 


H „r = ^u a o- z + H 


g 2 11 1 

+ 2 (A + X 1 °" 


[a + a 1 j . 


(37) 


This expression is analogous to the RWA dispersive 
Hamiltonian in Eq.(|33]l, with an extra contribution of 
X in the coupling term. However, this Hamiltonian is 
not diagonal in the eigenbasis of H 0 due to a f2 and 
a 2 . Nevertheless, for g/co r < 1, we can reinterpret 
the result as the state-dependent shift of the resonator 
frequency’s potential curvature u 2 m. Hence, the 
dispersive Hamiltonian with non-RWA gives rise to a 
shift in the oscillator frequency of 

u r —> u r , nr - u r ± g 2 j— + — j , (38) 

which implies that dispersive readout is also possible 
even in the ultrastrong coupling regime. Looking 
back at Eq.(|30|), we notice that both the time-averaged 
coupling strength g ar and g ar contribute to the 
two-photons terms and a 2 in the Hamiltonian. 
The time average coupling strength associated with 
counter-rotating terms also leads to an extra qubit 
dependent shift g 2 cr z /A. 

Next, we study dynamics of the dispersive case 
from strong coupling regime to ultrastrong coupling 
regime use the time-averaged coupling strength in 
IRWA by numerical simulation. In Fig. [3] we show 
the frequency shift of the resonator as a function of 
normalized coupling strength g/u r for positive detuning 
A > 0 and negative detuning A < 0 in RWA (blue 
dashed-dotted lines), non-RWA (red dashed lines) and 
IRWA (black solid lines) with Gaussian weighting 
function. It is clear that the RWA results have a 
totally different trend compared with the non-RWA 
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Figure 3: (Color online) Resonator frequency shift with single qubit 
in dispersive case as a function of coupling strength g/o r obtained 
with RWA (blue dashed-dotted lines), non-RWA (red dashed lines) 
and IRWA (black solid lines) for (a) positive detuning A > 0 (A = lOg) 
and (b) negative detuning A < 0 (A = -lOg) with Gaussian weighting 
function of width cdk = 10|A|. The breakdown of RWA is obvious, 
where it underestimates the dispersive shift for positive detuning and 
gives rise to a spurious shift in the absence of qubit ( oj a = 0) for 
negative detuning. 


results, especially for larger coupling strength. It 
underestimates the dispersive shift for positive detuning 
and predicts a shift even when the qubit’s frequency co a 
tends to be zero for negative detuning (A —» co r as g 
increases). This indicates the breakdown of RWA in 
predicting the dispersive resonator frequency shift in 
ultrastrong coupling regime. Meanwhile, our IRWA 
shows the manifestation of the counter-rotating terms 
as the coupling strength increases. 

In Fig. [4] we show the frequency shift of the 
resonator as a function of detuning A/co r in RWA (blue 
dashed-dotted lines), non-RWA (red dashed lines) and 
IRWA (black solid lines) with coupling strength of 
g/co r = 0.1 and Gaussian weighting function of width 
(Ok - 10|A|. For this relatively large coupling strength, 
it is shown that the RWA results underestimate the 
resonator frequency shift for JCM, whereas the IRWA 
predictions agree quite well with the non-RWA results 
for QRM. 
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Figure 4: (Color online) Resonator frequency shift with single qubit 
in dispersive case as a function of detuning A/cj r obtained with RWA 
(blue dashed-dotted lines), non-RWA (red dashed lines) and IRWA 
(black solid lines) for (a) positive detuning A > 0 and (b) negative 
A < 0 with Gaussian weighting function of width ldk = 10|A| and 
gl(O r = 0.1. 



4.2. Dispersive Regime with multi-qubit 

We now extend our discussion to multiple qubits 
coupled to a single mode resonator, where the 
time-averaged Hamiltonian takes the form lf40l 

H 7abi = T X + fwJra ^ a + H Yj ( 8 ' X + + S J arY J +), 

j j 

(39) 

with X J ± = acr J + ± a^cr J _ and Y ] ± = acrL ± a^cr ] + . Applying 
the unitary transformation 


U nq = exp 


^(AjXi + AjYi) 


v J 


(40) 


where the last four terms are the effective coupling 
between the qubits mediated by the resonator. 

To illustrate, we now take a two-qubit system as 
an example. In RWA, the counter-rotating terms are 
discarded because gj7 ~ 0. By setting all g{ = g and 
u> J a - to a , we obtain 



and expanding the transformed Hamiltonian to the where the interqubit interaction is of isotropic XY type, 

second order in Aj and A ; , we obtain the dispersive <r J <r k ¥ + <r\ c^ , '. In a frame rotating at the qubit’s 
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frequency, H 2q ' d generates the evolution 


U 2q ’ d 


=exp 


-iJ r t\a^a + 




'10 0 Of 

0 cos 7,-f i sin 7,7 0 
0 i sin 7,7 cos 7,7 0 

,0 0 0 1 , 


(43) 


with /, being the identity operator in resonator space and 
the effective coupling strength being 

g~ 

Jr = K ~. (44) 

This has been employed to generate qubit-qubit 
entanglement and quantum gate operations j6][40]. For 
instance, by turning on the coupling for a period t = 
nA/4g 2 , we can generate a V/SWAP gate which can 
be used to transform the state |ey,g,t) into an entangled 

state 1/ ^2(\ej,g k ) + i\gj,e k )). Here, \ef) and | g k ) are 
the excited state for y-th qubit and ground state for the 
k-th qubit, respectively. 

In ultrastrong coupling regime without RWA, all the 
terms will be retained. By setting all g J r — g J ar = g and 
a>i = cj a to be equal, we obtain 



where the interqubit interaction is of Ising type cr J x cr k x . 
The evolution operator reads 


t r2 q,d 

U n r = exp 


-iJnrflt ^j (o- J z + of) 


(cos J nr ^t 0 0 i sin J nr jf 

0 cos J m ,\t i sin 7„,jf 0 

0 i sin 7„ r jf cos7„ n if 0 
j sin J nr 1 1 0 0 cos J nr \t, 


'In 


(46) 


in the frame rotating at the qubit’s frequency with the 
effective coupling strength being 


Jnr, 0 — g 

Jnr,\ = g 2 



(47) 

(48) 



Figure 5: (Color online) The effective coupling strength in two-qubit 
dispersive case as a function of coupling strength g/co r obtained with 
RWA (blue dotted lines), non-RWA (red dashed lines) and IRWA 
(black solid lines and green dashed-dotted lines) for (a) positive 
detuning A > 0 (A = lOg) and (b) negative detuning A < 0 
(A = -lOg) with Gaussian weighting function of width cjk = 10|A|. 


It is worth noting that the extension from Eq.(|42|) to 
Eq.(|45| is not just a renormalization of the parameters. 
The effective qubit-qubit interaction type is indeed 
different, which will be clearer when we compare 
the evolution operators for RWA and non-RWA 
(Eq.(43i and Eq.(|46|)), where one is isotropic XY 
interaction while the other is Ising type interaction 
respectively. To understand this apparent sudden 
transition between RWA and non-RWA in dispersive 
regime for the multi-qubit, we invoke the time-averaged 
IRWA interpretation. From the effective Hamiltonian in 
Eq.gT}. the evolution operator can be written as 


u i =exp 


-iJirfit\a^a + - j (crl + of) 


< cos Ji r gt 0 0 i sin Ji r gf 

0 cos 7,>, it /sin 7;,-if 0 

0 / sin J ir \t cos 7,-,. p 0 

\.i sin Jirgt 0 0 COS Jirgt, 


Hr, 

(49) 


where the effective coupling strengths being 


j _ ( gi) (girf 

Jir, 0 - —;- 1 ---, 


Jir, 1 = gigr I “T + I ~ Sarglr I ^ ^ 


Jir,2 = grgar I ~ “ W I + Strgr | ^ ~ | ■ 


J n k 


J 


(50) 

(51) 

(52) 


In Fig. 0 we show the transition of the difference 
for the qubit-qubit interaction type from RWA case to 
non-RWA using the IRWA with Gaussian weighting 
function of width u>k = 10|A|. As we can see, when 
the coupling strength is very small, glu) r ^ 0.1, the 
IRWA curves are closer to the RWA curves for 7,>,i, 
whereas they are essentially zero for J ir g. When the 
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coupling strength is increased, the counter-rotating term 
coupling strength g ar starts to become significant, and 
hence, leads to the correction of /,,-j and manifestation 
of J ir 2 - As the coupling strength further increases, our 
IRWA curves start to deviate from the RWA curves and 
agree better with the non-RWA curves, showing the 
transitions from RWA to non-RWA. Eventually as the 
coupling strength reaches ultrastrong coupling regime, 
we regain the non-RWA results as in Eq. ((46]). 


5. Conclusion 


In this paper, we introduced the IRWA that is based 
on the time-averaging method for better understanding 
of the roles of the “counter-rotating” terms in the 
QRM and the transition between strong coupling and 
ultrastrong coupling regimes. The eigenenergies of 
the system were studied by combining the perturbation 
theory and IRWA for near resonance case. The results 
agreed well with the JCM predictions for small coupling 
strength, i.e. glu> r up to 0.1 and with the QRM 
results for larger coupling strength, i.e. g/o),- up 
to 0.3. We also showed that in dispersive regime, 
our IRWA predication gave a good explanation of the 
qubit-dependent frequency shifts in the single qubit 
scenario. This approach revealed the emergence of 
counter-rotating terms in the interqubit coupling, which 
leads to both quantitative and qualitative differences in 
the interaction strength and interaction type. Compared 
with other approaches liTTl Il2l 16] 171, our IRWA 
method allows us to gain the physical consequences 
of the co-rotating and counter-rotating coupling terms 
individually by tracing those terms separately. As 
a remark, there are several aspects that still can be 
explored with the idea of IRWA. For instance, by 
relating the measurement interval to the width of the 
time-averaging function u> K in our analysis, we can 
extend the result in ED to observe the transition from 
quantum Zeno effect to quantum anti-Zeno effect Il42l 
|43l . Our IRWA approach could also be applied to the 
studies of applicability of RWA in various phenomena, 
such as Berry phase in quantum systems El SI, 
asymmetric couplings IfiZl . and generalized multi-qubit 
quantum Rabi model lf32l . The IRWA might be useful 
as well in studying the dynamics of multiple coupling 
regimes in a single system, for example, by having one 
qubit coupled strongly in RWA regime and the other one 
operated ultrastrongly beyond RWA regime. 
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